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Abstract. The straylight contamination due to the 
Galactic emission (GSC, Galaxy Straylight Contamina- 
tion) entering at large angles from the antenna centre 
direction may be one of the most critical sources of sys- 
tematic effects in observations of the cosmic microwave 
background (CMB) anisotropies by future satellite mis- 
sions like Planck and MAP. For the Low Frequency In- 
strument (LFI) , this effect is expected to be particularly 
crucial at the lowest frequency channels. We describe 
here different methods to evaluate the impact of this ef- 
fect and compare it with other systematics of instrumen- 
tal and astrophysical origin. The results are presented 
in terms of simulated data streams and maps, Fourier 
series decomposition and angular power spectrum. The 
contributions within few degrees from the beam centre 
dominate the GSC near the Galaxy plane. The antenna 
sidelobes at intermediate and large angles from the beam 
centre dominate the GSC at medium and high Galactic 
latitudes. We find a GSC peak at ~ 15/iK and a GSC 
angular power spectrum above that of the white noise 
for multipoles I ;$ 10 albeit smaller than that of CMB 
anisotropies by a factor larger than ~ 10. At large mul- 
tipoles, the GSC affects the determination of CMB an- 
gular power spectrum significantly less than other kinds 
of instrumental systematics, like main beam distortions 
and 1// noise. Although the GSC is largest at low Galac- 
tic latitudes, the contamination produced by far pattern 
features at medium and high Galactic latitudes, peak- 
ing at ~ 4/iK, has to be carefully investigated, because 
the combination of low amplitude of Galaxy emission in 
those regions with the extremely good nominal Planck 
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sensitivity naturally makes high Galactic latitude ar- 
eas the targets for unprecedentedly precise estimation 
of cosmological CMB anisotropy. This paper is based on 
Planck LFI activities. 

Key words: Cosmology: cosmic microwave back- 
ground - Galaxy: general - Space vehicles - Telescopes - 
Methods: data analysis 



1. Introduction 

After the great success of COBE-DMR (Smoot 
et al. 1992, Bennet et al. 1996a, Gorski et al. 1996) 
which probed the gravitational instability scenario for 
structure formation through the detection of CMB 
anisotropies at few degree scales, and the recent balloon- 
borne experiments at high sensitivity and resolution on 
limited sky regions (De Bernardis et al. 2000, Hanany 
et al. 2000), supporting a universe model with Q to t ~ 1 
(Lange et al. 2000, Balbi et al. 2000, Jaffe et al. 2000), ul- 
timately, future progresses of the CMB anisotropy cos- 
mology will be based on two space missions, MAP (Mi- 
crowave Anisotropy Probe) (see Bennet et al. 1996b) by 
NASA and Planck by ESA, planned to be launched 
respectively in the years 2001 and 2007. 

In particular, the Low Frequency Instrument (LFI, 
Mandolesi et al. 1998) and the High Frequency Instru- 
ment (HFI, Puget et al. 1998) on-board Planck will 
cover together a wide frequency range (30-^900 GHz) 
which should significantly improve the accuracy of the 
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subtraction of foreground contamination from the pri- 
mordial CMB anisotropy, providing at the same time 
a gold mine of cosmological as well as astrophysical 
information (e.g. De Zotti et al. 1999a and references 
therein). 

To fully reach these scientific goals, great attention 
has to be devoted to properly reduce and/or subtract all 
the possible systematic effects. Detailed simulation codes 
have been developed and are continuously implemented 
to analyse the impact of several classes of instrumen- 
tal effects related to the behaviour of the optics, instru- 
ments and environment for a wide set of possible scan- 
ning strategies (e.g. Burigana et al. 1998, Delabrouille 
1998, Maino et al. 1999, Mandolesi et al. 2000a). Ulti- 
mately this effort leads to the optimization of the mis- 
sion design and the production of realistic data streams 
and simulated maps for data analysis tools as well as a 
realistic definition of Planck scientific performance. 

In particular, the behaviour of the Planck antenna 
patterns both at intermediate and large angles from the 
directions of beam centres have to be carefully consid- 
ered. The requirement on the rejection of unwanted ra- 
diation coming from directions far from the optical axis 
(straylight) is stringent for Planck and does not pertain 
only the telescope itself, but the entire optical system, 
including also solar panels, shielding, thermal stability 
and focal assembly components. The primary sources of 
error for the LFI are those due to imperfect off-axis re- 
jection by the optical system of radiation from the Sun, 
Earth, Moon, planets, Galaxy and the spacecraft itself 
(de Maagt et al. 1998). The variations of the spurious 
straylight signal during the mission introduce contam- 
inations in the anisotropy measurements. 

The antenna response features at large angular scales 
from the beam centre (far sidelobes) are determined 
largely by diffraction and scattering from the edges of the 
mirrors and from nearby supporting structures. There- 
fore they can be reduced by decreasing the illumination 
at the edge of the primary, i.e. increasing the edge taper, 
defined as the ratio of the power per unit area incident 
on the centre of the mirror to that incident on the edge. 
Of course, the higher is edge taper, the lower is the side- 
lobe level and the straylight contamination. On the other 
hand, increasing the edge taper has a negative impact on 
the angular resolution for the fixed size of the primary 
mirror (e.g. Mandolesi et al. 2000a). A trade off between 
angular resolution and straylight contamination has to 
be found. 

The main aim of this work is to evaluate the im- 
pact of the Galactic emission as a source of straylight 
for Planck LFI. We will then compare it with the ef- 
fects generated by other kinds of systematics, the main 
beam distortion introduced by optical aberrations and 
the 1// noise related to gain fluctuations in LFI radiome- 
ters, and with the astrophysical contamination from the 



Galaxy and the extragalactic sources in the main beam. 
At LFI frequencies, the Galaxy straylight contamina- 
tion (GSC) is expected to be particularly crucial at the 
lowest frequencies, due to the increase of synchrotron 
absolute emission and anisotropies with the wavelength. 
For simplicity, we limit our analysis to the case of the 
30 GHz channel, but the methods presented here can be 
extended to higher frequencies. 

In Sect. ^ we briefly describe the basic recipes of our 
simulation code, discussing the geometrical aspects rele- 
vant for the full sky convolution, the format of the opti- 
cal simulation output performed by the ESA (de Maagt 
et al. 1998) that are adopted in the present work, the 
conversion from data streams to maps and the compu- 
tation of the Fourier modes and of the angular power 
spectra, two different estimators of the GSC impact. We 
estimated the expected GSC on the basis of the antenna 
integrated response from angular regions at different an- 
gles from the beam centre, and on the level of Galaxy 
emission. In Sect. [| we focus on the integration accu- 
racy of our computations and test the consistency of the 
code by assuming simplified input maps and antenna 
patterns. The main results concerning the evaluation of 
the GSC are presented in Sect. ||. In Sect. || they are 
compared with the effects introduced by other kinds of 
instrumental effects and several sources of astrophysical 
contamination. Finally, we draw our main conclusions in 
Sect. §. 

2. Simulation of Planck observations 

The selected orbit for Planck is a Lissajous orbit 
around the Lagrangian point L2 of the Sun-Earth sys- 
tem (e.g. Mandolesi et al. 1998). The spacecraft spins at 
1 r.p.m. and the field of view of the two instruments 
is 5° around the telescope optical axis at a given angle 
a from the spin-axis direction, given by a unit vector, 
s, choosen outward the Sun direction. In this work we 
consider values of a ~ 80° — 90° 111. The spin axis will 
be kept parallel to the Sun-spacecraft direction and re- 
pointcd by ~ 2.5' every ~ 1 hour. In addition, a pre- 
cession of the spin-axis with a period, P, of ~ 6 months 
at a given angle /3 ~ 10° about an axis, /, parallel to 
the Sun-spacecraft direction (and outward the Sun) and 
shifted of ~ 2.5' every ~ 1 hour, may be included in 
the effective scanning strategy. These kinds of scanning 
strategies do not modify the angle between the spin axis 
and the Sun-spacecraft direction, avoiding possible ther- 
mal fluctuations induced by modulations of the Sun il- 
lumination, and allow to achieve nearly full or full sky 
coverage. Hence Planck will trace large circles in the 
sky. The detailed distribution in the sky of the num- 
ber of observations per pixel depends on the adopted 
scanning strategy, the telescope design and the arrange- 

1 Most recently it was recommended to choose a = 85°. 
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ment of the feed array on the telescope focal surface. 
The scanning strategy and spacecraft geometry have to 
be carefully addressed in order to minimize systematic 
effects before and after the data analysis and for ensur- 
ing that the sky coverage be as complete and spatially 
smooth as possible. 

The code we have implemented for simulating 
Planck observations for a wide set of scanning strate- 
gies is described in detail in Burigana et al. (1997, 1998) 
and in Maino et al. (1999). Here we consider simple scan- 
ning strategies, namely with the spin axis kept always in 
the ecliptic plane. Under this assumption, the geometri- 
cal input parameters relevant for the scanning strategy 
are the angle a between the spin axis and the telescope 
line of sight and the beam location in the telescope field 
of view. 

In this study we neglect the small effects introduced on 
the GSC evaluation by the Planck orbit by simply as- 
suming Planck located in L2 and consider spin-axis 
shifts of 2°every two days and 180 samplings per scan 
circle. These simplifications allow to speed up the simu- 
lation without significantly affecting our understanding 
of the main effects introduced by the GSC, because of 
the decreasing of the Galaxy fluctuation level at small 
scales and because the effects of pattern features we want 
to study here occur at ~ degree or larger scales. In fact, 
a recent study of synchrotron emission (Baccigalupi et 
al. 2000a) based on recent high resolution surveys at 
low at medium latitudes indicates a steeper power law 
of the total intensity angular power spectrum where dif- 
fuse emission dominates. 

2.1. Coordinate transformations and input antenna pat- 
tern 

For the simple scanning strategies considered here, the 
vector s, which gives the spin axis direction, and a unit 
vector p, parallel to the telescope line of sight z, are sim- 
ple functions of the time and of three scanning strategy 
parameters: the spinning frequency f s , the angle a and 
the orientation of the Planck telescope line of sight at 
the beginning of each scan circle that we arbitrarily fix 
close as much as possible to the positive z ecliptic axis. In 
additon, we set the rotation of p clockwise with respect 
to the positive direction of s. 

On the plane tangent to the celestial sphere in the 
central direction of the field of view, i.e. on the field of 
view plane of the Planck telescope, we choose two co- 
ordinates x and y, with unit vector u and v respectively, 
according to the convention that the unit vector u points 
always toward s and that x, y, z is a standard cartesian 
frame, referred here as "telescope frame" □. 

2 We note that the "telescope frame" defined above may be 
or not equivalent to the "primary mirror frame" depending 
on the absence or presence of a tilt angle between the pri- 



In general, the beam centre will be identified by its 
unit vector b in the frame x, y, z or equivalently by the 
coordinates, 2:0, yen of its projection on the plane x 7 y 
or, more usually, by its corresponding standard polar 
coordinates, the colatitude 9 b and the longitude 4>b- The 
HFI feedhorns will be located in a circular area at the 
centre of the focal plane, and LFI feedhorns in a ring 
around HFI. Therefore the corresponding positions of 
LFI beams on the sky field of view are significantly off- 
axis. For a telescope with ~ 1.5 m aperture the typical 
100 GHz LFI beam is located at ~ 2.8° from the optical 
axis, whereas the 30 GHz beams are at about ~ 5° from 
it. 

The shape of the main beam computed by ESA (de 
Maagt et al. 1998) is provided in a regular equispaced 
grid on x, y about the beam centre. We can then perform 
the convolution of the main beam (within a choosen an- 
gle from the beam centre) with the sky signal directly in 
this frame. 

According to the standard output of the "GRASP8" 
code for optical simulations as performed de Maagt et al. 
(1998), we describe the antenna pattern response, J, at 
large angles from the beam centre by using two standard 
polar coordinates 9bf (between 0° and 180°) and <fibf (be- 
tween 0° and 360°) referred to the "beam frame". This 
corresponds to the standard cartesian "beam frame" 
Xbf,ybf,Zbf which is obtained by the "telescope frame" 
x, y, z when the unit vector of the axis z is rotated by an 
angle 9b on the plane defined by the unit, vector of the 
axis z and the unit vector b up to reach b U. We use here 
the polar coordinates 8bf,4>bf for the convolution of the 
antenna pattern with the sky signal at large angles from 
the beam centre. For the antenna pattern at "interme- 
diate" (namely up to few degrees from the beam centre) 
and "far" (namely for the entire solid angle) angular dis- 
tances from the beam centre different equispaced grids 
in 9b f are available: more refined for the former, because 
the response variations are stronger close to the beam 
centre, and less refined for the latter, where the relevant 
response variations occur on degree or larger angular 
scales. 

The orientation of these frames as the satellite moves 
is implemented in the code. For each integration time, 
we determine the orientations in the sky of the telescope 
frame and of the beam frame and compute the pattern 



mary mirror axis and the central direction of the sky field of 
view. For example, the current optical configuration forseen 
for Planck exibits a tilt angle of ~ —3.75°. 
3 Properly, the standard of "GRASP8" code for optical sim- 
ulations is to provide the far pattern in terms of a "colati- 
tude" angle 0' bf (between -180° and 180°) from the "polar" 
axis parallel Zbf and an "azimuthal" angle (j>'bf (between 0° 
and 180°) related to 9 b f and cftbf by: <j)'bf — 4>bf if <fo/ < 180° 
and <j}' bf =<f>bf- 180° if <j> hf > 180°; 0' bf = 6 bf if <j> bf < 180° 
and e' bf = -6 b f if <t> bf > 180°. 
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response in each considered sky direction, thus perform- 
ing a direct convolution with the sky signal. 



2.2. Optical results 

The telescope design of the Carrier Configuration for 
Planck - the current baseline in which Planck and 
FIRST are lauched together and separate in orbit - is 
based on the Phase A study. The configuration is a off- 
axis Gregorian Dragone-Mizuguchi telescope with the 
main reflector oversized at 1.5 meter projected aper- 
ture for reducing the spillover at the primary mirror. 
The subreflector axis is tilted at 14° with respect to the 
main reflector axis. The off-axis design of Planck tele- 
scope (see Fig. |l|) introduces particular features in the 
full sky antenna response. In Fig.|l| a sketch on the sym- 
metry plane of the telescope is reported. For simplicity, 
a fccdhorn is located at the centre of the focal plane. Its 
pattern on the sky is along the direction of the optical 
axis of the telescope (see bottom panel of Fig. |l|). The 
unwanted ray directions are shown in the upper panel of 
Figj^. The rays labelled with 1 are coming from the sky 
directly into the feed. The rays labelled with 2 are those 
scattered into the feed by the subreflector only. These 
identified regions are the spillover past the subreflector 
for rays 1 and the spillover past the main reflector for 
rays 2. Considering the Planck Carrier Configuration 
(de Maagt et al. 1998), the rays lb and 2b are blocked by 
the spacecraft shields, which redirect the rays in an an- 
gular region close to the main beam. In Fig. || we show 
a map of the considered full pattern; see the connec- 
tion beetween the regions marked in Fig. |l| and Fig. ^. 
In Fig. H we show also several cuts at constant azimuths 
<fibf from 0°to 360° (from the bottom to the top). 

We have conservatively considered the worst case (de 
Maagt et al. 1998) for what concerns the straylight ef- 
fect: we use the antenna pattern computed at 30 GHz, 
the channel with the highest spillover and with the high- 
est Galaxy signal. We also included the shields for the 
Carrier Configuration. 

The main feature is the spillover (2a) at about 90° 
from the main beam (see also Fig. |^, where the main 
beam is located at North pole) which shows a response 
of w —60 dB with respect to the maximum extending for 
few tens of degrees in 9b f and in ipbf (around <f>bf = 0°, 
i.e. always quite close to the direction of the axis x in the 
"telescope frame"). Another relevant feature is the sub- 
reflector spillover (la), with similar response level and an 
angular extension close to the main beam (6b f ~ 30°) as 
we can see also in Fig.0. Other features are located on 
the northern semisphere also, due to the shields which 
block rays coming from the southern part of the sphere 
(sec de Maagt et al. 1998 for a more detailed discussion 
of their connection with the optical configuration) . 

The pattern has been calculated by using the Physi- 
cal Theory of Diffraction (PTD). Since pattern responses 
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Fig. 1. Sketch of the telescope design in the symmetry plane 
(<f)' b f ~ 0). This figure shows the carrier telescope configu- 
ration: a gregorian Dragone-Mizuguchi design with the pri- 
mary mirror enlarged to 1.5 meter for reducing the main 
spillover (2a). The upper panel shows the rays coming from 
unwanted directions and focussed into the feed: (la) spillover 
past the subreflector; (lb) spillover past the subreflector redi- 
rect by the shields; (2a) spillover past the main reflector; (2b) 
spillover past the main reflector and redirect by the shields. 
The bottom drawing shows the telescope with rays coming 
from the direction of the optical axis and focussed in a feed 
located at the centre of the focal plane. The spillover blocked 
by the shield, (lb) and (2b), is redirected close to the main 
beam as showed on the full-antenna pattern in Fig. |i| 

at levels smaller than about —60 dB are hard to measure, 
this is the most accurate method for predicting the side 
lobe response on the antenna. The validity of the simu- 
lations at very small levels of the sidelobes will be tested 
by measuring the antenna response of a fully represen- 
tative copy of the Planck telescope both in compact 
range and in outdoor far field test range facilities. 
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Fig. 2. Full antenna pattern, with response normalized to 
the maximum, for the carrier configuration. The color table 
is linear in dB and the true directivity at the maximum is 
49.36 dB. Pattern regions related to particular optical struc- 
tures are identified. See also the text. 



Fig. 3. Cuts of full antenna pattern in dB for the carrier 
configuration. The lines refer to the antenna response for 9b f 
between 0° and 180° and, from the bottom to the top, with 
4>bf a t steps of 30° from 0° to 360° (each cut is vertically 
shifted for graphic purposes to have the pattern response 
value at 9b f = 0° equal to the considered value of 4>bf)- 



2.3. Simple estimates 

Taking the level of Galactic emission and the antenna in- 
tegrated response from angular regions at different an- 
gles from the beam centre, we can provide first order 
estimates of the expected GSC. 

The region with dbf between ~ 1.2° and ~ 2° con- 
tains about 0.5% of the integrated response; the region 
between ~ 1.2° and ~ 5° contains about 0.6% of the 
integrated response and all the rest of the far pattern 
(0bf ~ 5°) contains about 1% of the integrated response. 
Of course, the remaining main integrated response falls 
in the "main" beam (up to ~ 1.2°) (see also Sect. 3.1 for 



a discussion on the choice of these characteristic angles). 
In addition, in the main spillover (2a) enters w 0.1-^0.2% 
of the integrated response. 

The sky signal at 30 GHz is known with pixel size 
of about 2.5° by COBE-DMR. For the present study at 
30 GHz the relevant astrophysical source is the Galaxy 
emission. We have implemented "small" angular extrap- 
olations (see e.g. Burigana et al. 2000a for further de- 



tails) for generating Galaxy maps with resolution of 
about 1.2°, corresponding to Quad-Cube resolution 7. 

For simple estimates, we note that in the adopted 
30 GHz Galaxy map there are ~ 13D with a signal (in 
terms of antenna temperature T a ) larger the 2 mK, ~ 
73D with T a > 1.5 mK and - 230D with T a > 1 mK, 
while the minimum signal is ~ 0.05 mK and about the 
50% of the sky shows a signal ~ 0.1 mK. 

By combining these numbers with the percentages 
of integrated responses falling within the above differ- 
ent angles from the beam centre, we expect to find a 
contamination peaking at about 10 /xK from the pat- 
tern regions between ~ 1.2° and ~ 5° and at few fiK 
from the pattern regions outside ~ 5°. In particular, in 
the main spillover (2a) we expect a signal peaking at 
~ 2 fiK when it looks at high signal Galactic regions. 
Similar contributions are expected from the pattern fea- 
tures at few tens of degrees from the beam centre. Of 
course, smaller contaminations (~ 0.5/aK) are expected 
when the relevant pattern features look at regions with 
low signal Galactic. 
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Numerical calculations, like those presented in the 
next sections, are required for more accurate estimates. 

2-4- From data streams to sky maps 

The input map is converted from its original Quad-Cube 
pixelisation to equal area, hierarchic HEALPix pixelisa- 
tion scheme (Gorski e t al. 1998), adopted in the present 
work (see also section |3.l| for details about the nominal 
resolution of this map). The final output of the simu- 
lation code relevant here are 2 matrices with a num- 
ber of rows equal to the considered number of spin- 
axis positions n s for one year of mission (n s = 180 
here) and a number of columns equal to the number 
of considered samplings along one scan circle (n p = 180 
here). In the first matrix, N, we store the pixel numbers 
corresponding to the main beam central directions for 
the considered 180x180 integrations; they are stored in 
HEALPix pixelisation scheme at n S id e = 32 (the num- 
ber of pixels ripi x in a full sky map is related to n s ide 
by n P i X — 12nf ide , Gorski et al. 1998). In the second ma- 
trix, G, we store the antenna temperatures "observed" 
by the considered portion of the antenna pattern for the 
above pointing positions. We neglect here the receiver 
noise and all the other systematics. These data streams 
are the first output of our simulations; they give imme- 
diately the impact of GSC and are useful to understand 
how this effect changes during the mission. 

From these data streams it is quite simple to obtain 
observed simulated maps, that can be visualized for ex- 
ample in mollweide projection: we make use of N and 
G to simply coadd the temperatures of those pixels ob- 
served several times during the mission. In this way we 
attribute to each pixel the average of the signals observed 
when the antenna pattern, due to the scanning strategy, 
is differently oriented in the sky and Galactic regions 
with very different signal intensities enter in the inter- 
mediate/far sidelobes; of course, by coadding different 
samples of the same location into pixels the systematic 
error per pixel is smaller than the systematic error in the 
most contaminated sample. This is because for different 
samples of the same location on the sky the sidelobes 
are pointing towards different regions of the sky, some 
brighter some fainter. 

2.5. Power spectra 

We can analyse both data streams and maps in terms of 
power at different scale-lengths or multipoles (or modes). 

In order to analyse separately each scan circle of sim- 
ulated data streams we follow the approach suggested by 
Puget & Delabrouille (1999) and decompose the time se- 
ries from the scan circle in Fourier series: 

m— n p /2 — 1 

3(0 = X! a m \/2cos(ra£ + £ ,m) , (1) 

m=0 



where g(£) is the signal at the position identified by an 
angle £ (between and 2ir) on the considered scan cir- 
cle. The 2 x n p coefficients a m and (j ra can be easily 
computed from the time series by solving a linear sys- 
tem. The amplitude, a m , of each mode m, analogous to 
the multipole I of a usual spherical harmonic expansion, 
has to be compared with that of the wanted signal and 
of other sources of noise. This approach is particularly 
interesting for the data analysis during the mission, as 
Planck will continuously scan different circles in the 
sky. 

For analysing the maps of coadded signals we use 
the standard approach of computing their angular power 
spectrum. We produce maps in the HEALPix pixelisa- 
tion scheme (Gorski et al. 1998) which takes advantage 
from the isolatitude of the pixels for a quick generation 
of a map from the coefficients ag m of the spherical har- 
monic expansion and vice versa (Muciaccia et al. 1997). 
[We will show the angular power spectra in terms of 
ST e (v) = ^£{2£+l)Ci(v)/4ir}. This is a very significant 
test for evaluating the GSC impact on Planck science, 
as the estimation of angular power spectrum of CMB 
fluctuations from the sky maps is one of the main ob- 
jectives of Planck mission. It offers also the possibility 
of directly comparing the GSC with astrophysical con- 
taminations and other sources of instrumental noise. Of 
course, it produces a "global" estimate of GSC effect, 
useful in the analysis of the whole sky maps. 

3. Accuracy and tests 

The antenna pattern is theoretically known from opti- 
cal simulation codes at the desired accuracy and res- 
olution compatible with the available computing time. 
The adopted grids have a resolution much better than 
those of currently available sky maps at Planck fre- 
quencies, although extrapolations both in frequency and 
in angular resolution of existing maps allow to produce 
more refined simulated maps (e.g. Burigana et al. 1998, 
2000a). On the other hand, the knowledge of Galactic 
emission at high resolution is not particularly relevant 
here and only small angular extrapolations up to a res- 
olution of about 1° has been implemented. Given the 
currently available input map resolution and the adopted 
simplified scanning strategies, we are able to derive the 
power of GSC only up to I ~ m ;$ 80, a mode/multipole 
range satisfactory for the study of the smooth features in 
the intermediate and far pattern response. A very good 
agreement is found with the high resolution computa- 
tions based on faster Fourier expansion methods (Wan- 
delt & Gorski 2000) which will be practically necessary 
for whole sky convolutions at angular resolutions higher 
than those adopted in this work. On the other hand, 
at larger multipoles the GSC effects are dominated by 
other systematics, as discussed in Sect. 5. 
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3.1. Sky, pattern grids and numerical accuracy 

We are interested in producing signal data streams as 
they would be observed separately by different angu- 
lar regions of the antenna pattern in order to under- 
stand the effect of the different pattern features as they 
project in the sky during the Planck observations. We 
have considered here three regions: 0° < 9bf < 
(main pattern); 9bf,i < &bf < (intermediate pat- 

tern); 9b}, 2 < 9b f < 180° (far pattern). We have choosen 
9bfs — 1-2° and 9^,1 = 5°. Of course, the choice of 9bf,i 
and #6/, 2 has to be appropriate to the considered antenna 
pattern: for a given telescope design it depends mainly 
on the considered frequency and only weakly on the ex- 
act feed location on the focal surface. For the 30 GHz 
channel, the main beam can be accurately measured in 
flight through planets (Mandolesi et al. 1998, Burigana 
et al. 2000b) up to response levels of « — 30 dB with re- 
spect to the peak response; 9b / ^ 1.2° corresponds to 
antenna responses lower than —40 dB, where the beam 
response probably becames highly difficult to measure in 
flight; 9bf\2 ~ 5° roughly divides pattern regions where 
significant response variations occur on angular scales 
less than 1° from those where they occur on ~ degree or 
much larger scales. 

The observed antenna temperature is given by 

_ fJ(9 bf ^bf)T a (9bf,<j> bf )dn 

i a , obs - jj( 0bf><kf)dSl ' w 

where J and T a are the antenna response and the sky 
antenna temperature in the direction given by 9b f , 06/ . 

The convolution of antenna response with the sky 
and the integration of the antenna pattern is simply 
computed by adding the contributions from all the pix- 
els within the considered solid angle, at resolutions cor- 
responding to n si de = 1024,512,64 respectively for the 
main, intermediate and far pattern in order to take accu- 
rately into account the pattern response variations. The 
main pattern is given in equispaced cartesian coordinates 
with Ax = Ay = 0.0005. The intermediate pattern and 
the far pattern are provided in equispaced "GRASP8" 
polar grids with A<j>' bf = 10° and A9' bf = 0.1° and 0.5°, 
respectively. When we extract all the pixels in the sky 
that contribute to the convolution within the consid- 
ered solid angle, the exact central position of each pixel 
typically does not coincide with a grid point where the 
pattern is known. Simple standard bilinear interpolation 
(e.g. Press et al. 1992) on the pattern grid has been im- 
plemented: this is fast, robust and accurate enough for 
the present purposes. 

An estimate of the error introduced by the above 
discretizations and interpolation/computation methods 
can be provided by comparing the convolutions obtained 
with different values of n si d e , for example by increasing 
n si de to 1024 for the intermediate pattern convolution. 



The numerical error is negligible ( ^ 1/iK, £ 0.5^K or 
0.2/zK for the convolutions with the main, intermedi- 
ate and far pattern, respectively). 

3. 2. Tests with schematic skies and patterns 

Checking the consistency of the part of the code that 
computes the signal entering in the main beam and in 
the intermediate pattern is quite direct: we simply ex- 
pect that the maps extracted from the corresponding 
data streams are respectively very similar or roughly pro- 
portional (according to the fractional signal entering at 
1.2° < 9 b f < 5°) to Galactic emission pattern, except 
for the beam smoothing. 

Testing the validity of the computation of the sig- 
nal entering the far pattern is not immediate, because 
it does not reflect in a simple way the Galactic emission 
pattern. We have verified the consistency of our simu- 
lation code by exploring simple cases for which we can 
easily predict the large scale symmetries of the maps de- 
rived from the data streams observed by the far pattern. 
We have assumed a simple antenna pattern, centred on 
the optical axis and perfectly symmetric in 4>b / , given by 
the sum of two gaussian shapes, one for the main beam 
and one for the main spillover located at 90° from the 
main beam centre, plus a constant low response level. We 
have performed tests with the following different very 
simple input skies: i) a spot at North Galactic pole: it 
produces a map with a well defined slab on the Galactic 
plane; ii) a slab on the ecliptic plane: the corresponding 
map shows a signal maximum at the ecliptic poles and 
decreasing toward the ecliptic plane; Hi) a slab on the 
Galactic plane: the corresponding map shows a signal 
maximum at the Galactic poles and decreasing toward 
the Galactic plane, where it exhibits small longitudinal 
modulations related to those of the solid angle subtended 
by the main spillover (this is due to the scanning sym- 
metry with respect to ecliptic coordinates and not with 
respect to Galactic ones). We have also verified that the 
angular power spectrum of the these maps presents a 
main peak at the multipole I = 2 and secondary peaks 
at its harmonic frequencies, as expected from the 90° 
symmetry of the adopted far pattern. 

4. Simulations results 

We have considered two options for the Planck Car- 
rier configuration. The first is exactly that considered 
by de Maagt et al. (1998) with a = 80° and including 
shields. In the second case we have used the same opti- 
cal results but with a = 90°: the corresponding results 
are then only indicative, being not perfectly consistent 
because the spacecraft design would be slightly differ- 
ent for this configuration; on the other hand, this case 
is instructive because it allows to start addressing the 
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question of the dependence of GSC on a basic parame- 
ter of the scanning strategy. We have considered the an- 
tenna pattern at 30 GHz with the beam centre located 
at 9 B = 5.62° and (f> B = 126.03°. 

4-.1. Analysis of the scan circle data streams 

Figures [| and || show the absolute signals entering the 
main, intermediate and far pattern and their ratios for 
the data streams of three representative scan circles for 
a = 80° and 90°, respectively. 

Note that the signal entering the intermediate pat- 
tern is roughly proportional to that in the main beam: 
two main relative maxima typically appear, related to 
the two crossings of the Galactic plane. The signal from 
the far sidelobes exhibits a clearly different and shifted 
angular behaviour, although two main relative maxima 
are again typically present. These are mainly due to the 
contributions from the pattern features (la) in the cases 
of the left and central panels, and to the main spillover 
(2a) in the case of the right panels, as they cross the 
Galactic plane. Note in fact the displacement between 
the maximum signal from the main and far pattern, of 
few tens of degrees (several tens of degrees) for the left 
and central panel (right panel). 

We have applied the Fourier series decomposition 
(see Fig. ||) described by eq. (1) to the sum of interme- 
diate and far pattern (i.e. for #&/ > 1.2°) data streams 
from the scan circles shown in Figures | and |. 

The same decomposition has been applied to white 
noise data streams, computed according to the LFI sen- 
sitivity at 30 GHz (e.g. Maino et al. 1999) averaged over 
a number of scan circles that spans an ecliptic longitude 
length arc equal to the FWHM = 33', i.e. essentially 
the sensitivity corresponding to half year mission. The 
white noise power is above that of the signal entering at 
9bf > 1-2°, for practically all the modes m ^ 3, becam- 
ing « 10 times larger at m w 10; this is essentially due 
to the strong decreasing of Galaxy fluctuations at small 
angular scales. 

No significant differences are found by varying a from 
80° to 90°; for this reason and for sake of simplicity, we 
will show in what follows only the results for a = 80°, 
the angle for which the optical simulations have been 
appropriately performed. 

All simulated data streams for a 1 yr mission are 
shown in Fig. [?] for the case a = 80° (similar patterns are 
obtained in the case a = 90°). In the right panel, note 
the vertical high signal line at A ~ 270°, corresponding 
to the main spillover (2a), and the two high signal fea- 
tures, close to this line, at about (A,/3) ~ (180°, 0°) and 
(A, j3) ~ (360°, —160°), corresponding to the pattern fea- 
tures at few tens of degrees from the beam centre. Note 
also how the azimuthally asymmetric far pattern reflects 
in the large difference between the two halfs (along A 
axis) of the right panel of Fig. corresponding to the 
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Fig. 8. Maps of GSC in Galactic coordinates from the 
coadding of the simulated data of the whole year for the case 
q = 80°. [For graphic purposes, the pixels close to the ecliptic 
poles, non observed because of the adopted scanning strat- 
egy parameters, have been filled with values close to those of 
the adiacient pixels]. In the bottom panel, note the two arcs 
corresponding to the main spillover (2a) and the central con- 
taminations corresponding to the pattern features at few tens 
of degrees from the beam centre. They clearly correspond to 
the features in the right panel of Fig. 7. Note also the wide 
"stripe" connecting the two ecliptic poles, which corresponds 
to the sky region observed a single time in the current 1 yr 
simulation, owing to the considered effective scanning strat- 
egy- 



first and second six months of observation. This redun- 
dancy can be exploited for an efficient subtraction of 
GSC in the data analysis. 

4-2. Straylight contamination maps and angular power 
spectrum 



By coadding the data streams as described in Sect. 2.5 
we can obtain the corresponding maps. This is shown in 
Fig. 8 for the case a = 80°, by coadding the simulated 
data from the whole year. 

As apparent in Fig. 8, the map from intermediate 
pattern is roughly proportional to that derived from the 
main pattern; the relative intensities are roughly scaled 
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Fig. 4. Top panels: absolute signals entering the main (solid lines), intermediate (multiplied by 50 for graphic purposes, 
dotted-dashed lines) and far (multiplied by 100, dashed lines) pattern for the case a = 80° for three scan circles corresponding 
to spin axis longitudes of 0°, 120°and 270° from left to right. The angle along the scan circle, ranging between 0°and 360°, is 
set to at the beginning of the stream from each scan circle. Bottom panels: the ratios between the signals shown in the top 
panels, intermediate/main (solid lines), far/main (multiplied by 0.1, dotted-dashed line) and far/intermediate (multiplied 
by 5 • 10 -4 , dashed lines). 



by the fraction of integrated response entering the two 
portions of the antenna pattern. 

On the contrary, the sky "observed" by the far pat- 
tern is very different. The signal is higher close to the 
Galactic plane, because of the features in the antenna 
pattern within 10° — 20° from the main beam, and at 
about 90° from the Galactic plane, because of the signal 
entering the main spillover (2a). 

We have computed the angular power spectra of 
these GSC maps (see Fig. ||) and compared them with 
the theoretical angular power spectrum of the white 
noise for a single 30 GHz receiver and for four receivers 
and with a typical CMB anisotropy angular power spec- 
trum (a tilted - n p = 0.9 - power spectrum with stan- 
dard CDM cosmological parameters and approximately 
COBE normalized) and with that of Galaxy fluctua- 
tions as seen by the main pattern. The most important 
contamination in terms of angular power spectrum de- 
rives from the signal entering in the intermediate pat- 
tern when all the sky is considered; on the contrary, if 
we consider only the regions at |6| > 30°, more crucial 



for Planck main science, the largest contributions to 
the GSC power spectrum derive from the far sidelobes. 
In general, the GSC power spectrum is larger than the 
white noise one at low multipoles (t ^ 5) but their ra- 
tio becames less than ~ 1/10 at £ ~ 50 and decreases 
further at larger multipoles, due to their different depen- 
dence on I u. 

5. Comparison with other sources of noise 

Many other sources of contamination, both instrumental 
and astrophysical in origin, may affect Planck obser- 
vations. We approach here a first comparison among the 
effects introduced by some of these systematics. 



The increase of the power of GSC from far sidelobes at 
£ ^ 50, evident in the plot, is not an intrinsic effect of the 
GSC but is generated by the "stripe" corresponding to the 
sky region observed a single time (see also Fig. 8) in the cur- 
rent 1 yr simulation with a — 80° and off-axis beam. 
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Fig. 5. The same as in previous figure, but for the case a = 90°. 



5.1. Other kinds of instrumental noise 

The impact of main beam distortions introduced by 
optical aberrations on Planck measurements has 
been carefully studied in several works (e.g. Burigana 
et al. 1998, 2000a, Mandolesi et al. 1997, 2000a). Burig- 
ana et al. (1998) discussed the impact of the main beam 
distortions for the representative case of an elliptical 
main beam shape. In general, the absolute rms addi- 
tional noise, in the range of few /zK, increases with the 
beam ellipticity. The combined effect of main beam dis- 
tortions and of Galaxy emission fluctuations increases 
the additional error at ~ 30 GHz by a factor ~ 3 with re- 
spect to the case of the essentially pure CMB fluctuation 
sky at high Galactic latitudes, whereas it produces only 
a small additional effect at higher LFI frequencies. In 
addition, the combined effect of main beam distortions 
and extragalactic source fluctuations is found to be very 
small at all LFI frequencies (Burigana et al. 2000a) com- 
pared to the noise induced by beam distortions in the 
case of a pure CMB sky. Then, we focus here further on 
the impact of the main beam distortions on the determi- 
nation of angular power spectrum of CMB fluctuations 
by considering the idealized case of the above pure CMB 
fluctuation sky. The kind and the magnitude of optical 
distortions depend on the details of the optical design; 



for aplanatic configurations currently under study (e.g. 
Villa et al. 1998 and Mandolesi et al. 2000b) the typi- 
cal main beam shape is roughly elliptical owing to the 
strong reduction of the coma distortion. We computed a 
full year simplified simulation both for a pure symmetric 
gaussian beam with FWHM = a V81n2 = 30' and for an 
elliptical gaussian beam with axial ratio r — 1.3 and with 
the same effective resolution {^a x a y = a) of the sym- 
metric beam (r = 1). We shift the spin axis at steps of 5' 
and consider a step of 5' between two samplings on the 
same scan circle. We computed the difference between 
the maps obtained from the elliptical and the symmetric 
beam by coadding the corresponding data streams and 
calculate the angular power spectrum of this difference 
map in order to understand which range of multipolcs 
is mainly affected. As expected (see Fig. 10), this effect 
is particularly relevant at quite large multipoles, close 
to the CMB peak, where GSC power significantly de- 
creases: the magnitude, of course, is related to the value 
of r. 

From optical simulations we know that r typically 
increases with the distance from the beam centre. From 
the present simulations we infer that a value of r quite 
smaller than 1.3, say less than ~ 1.1 up to —3 dB from 
the centre and less than ~ 1.2 up to — 20 dB from the 
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Fig. 6. Fourier decomposition, see eq. (1), of the signal at 8tf > 1.2° from the scan circles of Fig. ^ (a = 80°, top panels) 
and ^ (a = 90°, bottom panels). From the left to right: spin axis longitudes of 0°, 120° and 270°. Dashed and solid horizontal 
lines represent the white noise for one and four receivers, respectively. 



centre, is good enough for avoiding significant contam- 
inations in the data, in agreement with the indications 
inferred on the basis of the approximations of Burigana 
et al. (1998) for the rms noise added by a main beam 
elliptical distortion. 

The 1/f noise due to amplifier noise temperature 
fluctuations induced by gain fluctuations in Planck LFI 
receiver and its dependence on the relevant instrumen- 
tal parameters has been studied by Seiffert et al. (1997). 
It introduces additional noise in Planck observations 
which shows as stripes in final maps (Janssen et al. 1996) 
owing to the particular Planck scanning strategy. 
We have recently carried out detailed studies (Maino 
et al. 1999 and reference therein) on its effect on Planck 
LFI measurements and on the efficiency of destriping 
algorithms based on the use of the crossings between 
different scan circles (Bersanelli et al. 1996, Delabrouille 
1998) for a wide set of Planck scanning strategies. We 
extend here their simulations by relaxing the hypothesis 
of symmetric beam to study the impact of main beam 
distortion into the destriping algorithm. In Fig. [ll] we 
show the angular power spectrum of the receiver noise 
before and after applying the destriping algorithm when 
we include also the above elliptical distortion, for a sim- 



ple scanning strategy with a = 90° and a beam location 
at 9b = 2.8°, 4>b = 45°, a "mean" choice regarding to 
the destriping efficiency (Maino et al. 1999). We find 
that the destriping efficiency is not significantly affected 
by the additional uncertainty introduced by the "sys- 
tematic" differences among the observed temperatures 
resulting from different orientations of the main beam 
at the crossing points of different scan circles. 

As evident by comparing Figs. 9 -j- 11, there is a 
crucial difference between the angular power spectra of 
GSC, main beam distortion induced noise and 1// noise. 
The GSC affects particularly the determination of CMB 
angular power spectrum at low multipoles, whereas main 
beam distortions are critical at large multipoles. The 1/f 
noise affect both high and low multipoles, but destrip- 
ing algorithms are particularly efficient in removing high 
multipole features in the power spectrum. It is clear that 
all these effects have to be reduced both via hardware 
and software. The 1/f noise can be reduced indepen- 
dently of the other two, its magnitude being related es- 
sentially to the instrument stability and to the scanning 
strategy (Maino et al. 1999). On the contrary, a com- 
promise has to be reached between GSC and main beam 
distortion noise, being both related mainly to the opti- 
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Fig. 7. Synthetic view of the data stream from all scan circles for a = 80° . The ecliptic coordinates A and f3 properly : 
here to the longitude and latitude of the telescope axis. [For graphic purposes, in this plot the range of /3 between —80° 
—240° refer to the second half of each scan circle]. Compare with Fig. 4. 



cal design. As discussed in Mandolesi et al. (2000a), for a 
given telescope design, their relative weight is controlled 
by the edgetaper. The optical design has to be optimized 
to find a trade off for reducing the combined impact of 
these two effects. 

5.2. Astrophysical contamination 

The impact of foregrounds on the primary cosmological 
goal of Planck mission has been extensively studied 
in literature for what concern both Galactic and extra- 
galactic contaminations, of discrete and diffuse origin; 
Planck itself will be a good opportunity for studying 
cluster physics, many classes of extragalactic and Galac- 
tic sources and the diffuse emission from the Galaxy 
(e.g. De Zotti 1999a and references therein). Many ap- 
proaches have been studied to separate the different 
components of the microwave sky and for deriving their 
angular power spectra (e.g. Tegmark & Esftathiou 1996, 
Hobson et al. 1998, Baccigalupi et al. 2000b, and refer- 
ences therein). 

We consider here the foreground impact on CMB sci- 
ence and their comparison with the effect of instrumen- 
tal systematics. In Figure fi^ we report several estimates 



of the angular power spectra of different astrophysical 
components. 

The Galaxy angular power spectrum is known to de- 
crease with multipole order £: we show here the power 
spectrum derived from the map observed by the adopted 
main pattern by cutting (lower thin dotted-dashed line, 
"CG') or not (upper thin dotted-dashed line, "WG") the 
region at |6| < 30° and the power spectra proposed by 
Tegmark & Esftathiou (1996) for free-free (thick dot- 
ted line, "FF"), synchrotron (upper thick dotted-dashed 
line, "S") and dust (upper thick dashed line, "D") emis- 
sion at relevant Galactic latitudes. We show also for com- 
parison the power spectra for synchrotron (lower thick 
dotted-dashed line, "s") and dust (lower thick dashed 
line, "d") as derived by Prunet et al. (1998) and Bouchet 
et al. (1998) for a sky patch at medium latitudes. Of 
course, Galaxy contamination strongly depends on the 
considered region. 

In terms of 5Tg, Poisson fluctuations from extra- 
galactic unresolved discrete sources increase approx- 
imately proportionally to the multipole I. We show 
here, separately for radiosource (thin solid lines, "R") 
and the far infrared galaxy (thin dashed lines, "IR"), 
the Poisson fluctuation power spectra predicted by 
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Fig. 9. Comparison between the angular power spectrum 
of GSC from different pattern regions and that of the CMB 
anisotropy (thick solid line) and of the receiver white noise 
(dotted lines; upper line: a single receiver; lower line: four 
receivers) . Galaxy fluctuation power spectrum as seen by the 
main pattern without map cuts (upper dotted-dashed line) 
and by considering only the regions at |6| > 30° (lower dot- 
ted-dashed line). Angular power spectra of GSC from the in- 
termediate pattern without Galactic cut (upper solid - green 

- line) and at |6| > 30° (lower dashed - green - line) and 
from the far pattern without Galactic cut (lower solid - red 

- line) and at \b\ > 30° (upper dashed - red - line). See also 
the text. 

Toffolatti et al. (1998) as recently revised by Toffo- 
latti et al. (1999), De Zotti et al. (1999b) and references 
therein, on the basis of current source counts and as- 
suming evolution models and spectra in agreement with 
current data, when sources above 1 Jy (upper curves) 
or 100 mJy (lower curves) are detected and subtracted. 
Of course radiosources dominate at low frequencies. We 
have taken into account here a gaussian (FWHM=33') 
beam smoothing in all cosmological and astrophysical 
angular power spectra and consequently neglected it in 
the receiver noise angular power spectrum. 

At low multipoles, Galaxy contamination is larger 
then instrumental effects, dominated by the GSC (and 
possible residual 1// noise); on the other hand, if the 
Galaxy emission and anisotropy can be modelled at few 



Fig. 10. Angular power spectrum introduced by an elliptical 
main beam distortion with r — 1.3 (lower - blue - line). See 
also the text. 

percent accuracy and at high Galactic latitudes, instru- 
mental effects can became comparable with the residual 
Galaxy contamination. 

If not appropriately taken into account in the data 
analysis, main beam distortions may introduce at high 
multipoles an additional contamination comparable with 
that introduced from radiosource fluctuations, after their 
subtraction at few hundreds mJy level. 

6. Conclusions 

We have studied the impact of GSC on Planck ob- 
servations at 30 GHz, by considering different and com- 
plementary evaluation approaches: absolute and relative 
quantification of the impact on scan circle data streams, 
Fourier decomposition of scan circle signals, computa- 
tion of maps of GSC and evaluation of their angular 
power spectra. These different methods allow us to fo- 
cus on different aspects of GSC. 

No relevant differences are found by varying the angle 
a between the spin axis and the telescope line of sight 
from 80° to 90°. 

Our simulations show that the GSC peaks at values 
of about 15//K (a value comparable with the sensitivity 
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Fig. 11. Angular power spectrum of the 1/f noise before 
(upper - red - line with blobs) and after (lower - green - line, 
close to the level of theoretical white noise power spectrum 
for a single receiver) the destriping for the case of an elliptical 
main beam. 

per pixel), mainly due to the signal entering at few de- 
grees from the beam centre. Such values are found in the 
regions quite close to the Galactic plane, where in any 
case the "direct" (i.e. observed by the main beam) con- 
tamination from the Galaxy prevents an accurate knowl- 
edge of CMB fluctuations, as it is next to impossible to 
remove the Galactic signals to accuracy better than ~ %. 
These large contamination values, although critical for 
CMB anisotropy measurements near the Galactic plane, 
are not crucial for the determination of Galaxy emission, 
which is several order of magnitude larger. 

By considering all the pixels in the sky, the typical 
values of GSC are less than the 50% of the white noise 
sensitivity per pixel. 

The most crucial contamination derives from the sig- 
nal entering in the far pattern, in spite of its peak values, 
of about 4/iK, smaller than those obtained for the inter- 
mediate pattern regions. In fact, although this effect does 
not seem to be very large in amplitude (indeed, being 
nominally subdominant to power spectrum of expected 
receiver noise) it does dominate the GSC at medium and 
high Galactic latitudes, which are critical regions for the 



Fig. 12. Angular power spectrum of different sources of as- 
trophysical contamination compared with the CMB angular 
power spectrum (thick solid line, labelled with "CMB") and 
that of the white noise power spectrum for four receivers 
(dotted line, "N"). See the text for the meaning of the other 
lines and labels. 



extraction of the best quality results on CMB anisotropy. 
This could be also critical for the Planck polarization 
measurements which will take advantage from the two 
patches close to the ecliptic poles where the sensitivity 
will be several times better than the average, according 
to the scanning strategy and the feed array arrangement. 

As expected on the basis of the behaviour of Galaxy 
emission angular power spectrum, the GSC affects the 
determination of the CMB angular power spectrum 
mainly in the low multipole region and much less at large 
multipoles, particularly when compared with the other 
instrumental effects considered here, the main beam dis- 
tortion and the 1/f noise. The additional noise intro- 
duced by the main beam distortion can be in principle 
subtracted in the data analysis, provided that the beam 
pattern is accurately reconstructed. 

Of course, a substantial improvement in the data 
analysis is necessary to jointly treat all the systemat- 
ics, of instrumental and astrophysical origin. From the 
telescope design point of view, the best optimization of 
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the edge taper requires a trade off between the main 
beam distortion and the GSC effects. 
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